Quantum state tomography of single photons in the spectral degree of freedom 
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Quantum State Tomography (QST) of optical states is typically performed in the photon number 
degree of freedom, a procedure which is well understood and has been experimentally demonstrated. 
However, optical states have other degrees of freedom than just photon number, such as the spatial 
and temporal/spectral ones. Full characterization of photonic states requires state reconstruction 
in these additional degrees of freedom. In this paper we present a technique for performing QST of 
single photon states in the spectral degree of freedom. This is of importance, for example, in quantum 
information processing applications, which typically impose strict requirements on the purity and 
distinguishabihty of independently produced single photons. The described technique allows for full 
reconstruction of the spectral density matrix, allowing the purity and distinguishabihty of different 
sources to be readily calculated. 

PACS numbers: 42.50.-p,42.50.Ar,42.50.Dv 



The preparation of single photon states is vital to the 
development of many optical quantum information pro- 
cessing protocols This has motivated much research 
into photon engineering techniques HHSSISHSIII 
In order to be useful for quantum information process- 
ing applications there are typically very stri ngen t re- 
quirements on the purity, distinguishabihty [1(1 HH , and 
mode-structure |l2( of prepared single photons. There- 
fore, characterizing single photon sources is of great prac- 
tical interest. 

The most general approach to characterizing quantum 
states is via Quantum State Tomography (QST) [l3j . 
which, using many identical copies of a given state, pre- 
scribes an approach for complete reconstruction of the 
density operator from the measurement of experimen- 
tally accessible observables. In the context of optical state 
characterization this is typicall y pe rformed in the photon 
number degree of freedom [lj, [lj, , most notably us- 
ing Optical Homodyne Tomography (OHT) [l^. I n this 
paper we describe an approach for performing QST of 
single photons in the spectral/temporal degree of free- 
dom. From the reconstructed density operator, important 
measures, such as the spectral purity and distinguisha- 
bihty of different sources, can be readily calculated. Pre- 
vious work has examined the issue of experimentally de- 
termining the temporal wave-packets of single photons 
|18| . However, our work is more general in that allows 
for complete reconstruction of the spectral density oper- 
ator. 

Conceptual overview — We begin with a concep- 
tual overview to provide some intuition into our pro- 
tocol. The main intuitive notion upon which our pro- 
tocol is based is that the relationship between a pho- 
ton's wavepacket in conjugate domains (e.g. time and 
frequency) provides information about the purity of the 



state in those degrees of freedom, and hence the off- 
diagonal elements of the density matrix. For example, 
in the limit of pure photons we approach 'transform lim- 
ited' wavepackets, meaning that the frequency and time 
domain wavepackets are directly related by Fourier trans- 
formation. In this case the state is of the form 

\ip) = / ip(u)a} (w) du\vac) = I ^(t)h\t)At\vac), (1) 



where a'(oj) (a T (i)) is the single frequency (time) pho- 
tonic creation operator, ip(t) is the inverse Fourier trans- 
form of ip(uj), and the integral run over all frequencies 
(times). For impure states this relationship breaks down. 
As an example, consider a photon which is a mixture 
of identical, but temporally displaced wavepackets. This 
might arise when a photon source exhibits 'time-jitter'. 
Such a photon can be described in the form 



P = 



/(kMtiWte)!*! - t c )(t 2 - t c \ dti df 2 dt c 



(2) 

where f(t) characterizes the mixture (i.e. time-jitter), 
ip(t) characterizes the temporal wave-packet of the com- 
ponents in the mixture, t c denotes a particular center 
time, and \t) = a^(t)\vac) is the time-specific single pho- 
ton state at time f. Measuring this state in the time do- 
main using time-resolving detectors we will observe the 
envelope of the entire mixture, roughly /(f), assuming 
ip(t) is narrow compared to /(f). However, measuring this 
state in the frequency domain with frequency-resolving 
detectors will not yield an envelope given by f(u>). This 
is because the Fourier transforms of each of the compo- 
nents in the mixture are identical up to irrelevant phase 
factors. Thus, in the spectral domain we observe the en- 
velope of the identical pure components of the mixture, 

i.e. ip(w). 

This observation provides a conceptual hint that in or- 
* rohde@physics.uq.edu.au URL: http : //www. physics .uq. edu. au/peoprk^tcbfiia/laasure the spectral density matrix we will need 
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to perform transformations in the conjugate domains, fol- 
lowed by a comparison of these transformed states. While 
this is a very qualitative description, our protocol is in 
fact a direct implementation of this idea - we create a su- 
perposition of an incident state across two paths, in one 
path implementing a frequency domain transformation, 
and in the other a time domain transformation, and fi- 
nally we implement a 'comparison' by interfering the two 
paths. 

Description of the protocol — We now describe 
the protocol in detail. We begin with a single photon 
state with completely arbitrary density operator in the 
spectral degree of freedom. This can be expressed in the 
form 



Pin 



p{u)\, 0J2 ) <v(u>i) I vac) {vac\a{u>2) dwi du>2 
p{u)\, U)2)\W\) {u)2\ dwi dcJ2- (3) 



Our goal is to experimentally sample arbitrary elements 
of the spectral density matrix p{ui\,uj2)- 

We employ a Mach-Zehnder interferometer and per- 
form different temporal and spectral operations in the 
two arms of the interferometer. The proposed layout is 
shown in Fig. ^ In the upper arm we induce a tunable 




FIG. 1: The proposed experimental setup. We employ a bal- 
anced Mach-Zehnder interferometer. In the upper arm is 
an Acousto-Optic- Modulator (AOM), followed by a phase 
shifter. In the lower arm is a time delay. We post-select upon 
measurements where a photon is detected at one of the out- 
puts. 

spectral shift using an Acousto-Optic-Modulator (AOM). 
The AOM implements the transformation 0] 



Ol.outM = V 1 - £«2M + iy^axiuj - 8) 
O2,out(w) = yl - £ai(w) + iy/£a 2 (u + 8), (4) 

where 8 is the modulation frequency applied to the AOM, 
and £ is related to the modulation frequency. d\ and 0,2 
denote distinct spatial modes. We assume the upper arm 
of the interferometer runs through mode 1 , while the vac- 
uum state is incident upon S2, and ci2,out is discarded. 
Following the AOM is a phase shifter, which implements 
the transformation 

o(w) -» e l6 a{uj). (5) 
In the lower arm is a temporal delay, implementing 

-» e-^ T b(uj). (6) 



The 50/50 beamsplitters are modeled by the transforma- 
tion 



flout M 
flout M 



a(u>) + b(u>) 



A/2, 



a(u) - b(u) /V2. 



(7) 



A simple calculation now shows that the conditional out- 
put state for detecting a photon in mode A can be ex- 
pressed in the form 



Pa(t, 8, 9) 



1 

4. 

+ e M - i ™ 1 p(wi>W2-*) 
+ e - ie + l ™V( Wl -o,w 2 ) 
+ p(u>i — 8, lu 2 — d)]\u)i)(u2\ dwj dcL>2- (8) 



The probability of detecting a photon at output mode A 
is just the normalization of this state. Thus, 



Pa(t,S, 



p(cj,cj- 6) 



-i6-\-iruj 



p{uj — 5, id) + p(uj — 5, id — S) du;, 

(9) 

where we have used the identity tr(|o;i)(a; 2 |) = (u>i\u>2) = 
5(u>i—u>2)' Now we can apply several identities to simplify 
this expression, 

p(oj, oj) duo — J p(oj — 6, lj — 8) dej = tr(pi n ) = 1, 

e^p(u; - S, w) duj = T-l T [p(u> - 8, cj)} , 

p(u> - 8,u) = p(iv,ui - 8)*, (10) 

where T~ x denotes an inverse Fourier transform and sub- 
scripts are used to explicitly denote the change of variable 
taking place. Thus, wc have 



Pa(t, 8,1 
P b (t, 8,1 



1 1 

2 + 2 

1 1 

2 ~ 2 



Rc[ e *^ T [p(c^-o)]] , 
Re[e ie T-i T [p(^-S)}}.(ll) 



Note that in this derivation we have assumed £ = 1, 
corresponding to no loss into the discarded mode in the 
AOM. When this is not the case this effectively unbal- 
ances the interferometer. This can easily be compensated 
for by introducing an equal loss into the other arm of the 
interferometer. While this reduces the overall count rate, 
it does not undermine the operation of the protocol since 
we can simply post-select upon events where a photon ar- 
rives at one of the output ports. Thus, in our calculations 
we are justified in restricting ourselves to the lossless case, 
and Eg. [TTl refers to the post-selected scenario. 

The procedure for reconstructing the input density 
matrix is as follows. For every setting of r and 8 we 
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perform two measurements, with 9 = and 9 = it/2, 
allowing both the real and imaginary components of 
•7"u)->"7-[p(k ; ) w — 5)] to be measured at every point. This 
function can then be inverse transformed to establish 
the diagonal cross section of the input density matrix, 
p{lo,lo — <5) V lo. By measuring such cross sections for 
different 5, the complete density matrix can be recon- 
structed. Several examples of spectral density matrices 
and their corresponding measurement statistics are illus- 
trated in Fig. 

Note that in the case where we do not introduce a fre- 
quency shift, <5 = 0, only the diagonal elements of the 
density matrix are accessible. This allows the spectral 
distribution to be measured, but not the coherences be- 
tween different spectral components. This is expected, 
and essentially corresponds to existing interferometric 
techniques for characterizing photon wave-packets. 

Tolerance against experimental imperfections 
— We now turn our attention to major sources of ex- 
perimental imperfection that are likely to arise in the 
implementation of this scheme. Firstly, our protocol is in- 
herently resilient against source and detector inefficiency. 
This is because the protocol is post-selected upon events 
where a photon is detected at one of the outputs. 

The second major source of error is likely to be spa- 
tial mode-mismatch at the second beamsplitter in the in- 
terferometer. We now consider the effect this has on the 
operation of the protocol. We model this by first expand- 
ing our representation for the single photon state to in- 
clude the transverse spatial degrees of freedom. We define 
new mode-creation operators A^Jui) and B^,Juj), 
where ip(x,y) is the transverse spatial wavefunction of 
the incident photon. These operators create photons at 
a particular frequency, but with arbitrary spatial distri- 
bution. Formally, 



4h 



i/j(x,y)a* l (x,y,uj)dxdy, (12) 



where aJ(x, y, ui) is the space- and frequency-specific cre- 
ation operator for mode A. Note that with this definition 
we assume the spatial distribution function to be fre- 
quency independent. We model mode-mismatch by ap- 
plying an arbitrary transformation to the spatial wave- 
function in one arm of the interferometer, A^(uj) — > 
A^i{lo). The calculation proceeds as before, and we em- 
ploy the identity 



A\p {uJi)\vac){vac\Aip,(uj 2 ) 



tr 



= 5(ui-u) 2 ) J J ip(x,y)*ip'(x,y)dxdy 

= 7<5(wi — w 2 ), (13) 

where 7 parameterizes the degree mode-mismatch (0 < 
7 < 1, with 7=1 corresponding to perfect spatial mode 
overlap, and 7 = to no mode overlap). It can now be 



shown that the output probabilities are given by 

Pa(t,6,9) = ~ + iRe [je ie ^-i T [p(^ ~ *)]] . 
Pb(t, 5,6) = \ - ^Re fre^i^, w - 6)}] .(14) 

Thus, spatial mode-mismatch simply reduces the visibil- 
ity of the probability fringes by a constant factor. This 
can be compensated for using some a priori knowledge 
of the behavior of P4 and Pg . When r = (5 = 6 = Owc 
expect all photons to exit through port A, as per an or- 
dinary balanced Mach-Zchndcr interferometer. Thus, by 
measuring the statistics of Pa(0, 0, 0) and Pb(0, 0, 0), the 
mode-mismatch parameter, 7, can be directly inferred. 
With knowledge of 7, the inversion procedure will faith- 
fully reproduce p- m . 

Variations of the scheme — There is nothing 
unique about our choice of transformations. Other com- 
binations of transformations in the interferometer arms 
that are known to work include spectral-shift/spectral- 
filter (upper/lower arms), temporal-filter/spectral-filter, 
and temporal-delay/temporal-filter. All of these combi- 
nations allow for complete characterization of the spec- 
tral density operator in principle. However, due to the 
difficulty in implementing sub-wavepacket filtering (espe- 
cially in the time domain), these variations are unlikely 
to be experimentally feasible. 

Experimental considerations — In order to suc- 
cessfully reconstruct the spectral density operator we 
need to be able to induce incremental spectral/temporal 
shifts over a range on the order of the spectral/ temporal 
bandwidth of the photon. While implementing temporal 
delays with such high resolution and over an essentially 
arbitrary scale is experimentally feasible, our ability to 
induce frequency shifts via AOM's is more limited. Cur- 
rent AOM's can induce frequency shifts on the order of 
GHz, implying that this technique is limited to photons 
with, at most, spectral bandwidths of this order. Unfor- 
tunately this rules out application of this technique to 
some widely used single photon engineering techniques, 
such as ultra- fast parametric down conversion, which in- 
herently have spectral bandwidths outside this range. 

Conclusion — We have described an approach for to- 
mographically reconstructing the density operator of sin- 
gle photon sources in the spectral degree of freedom. This 
could prove useful in fully characterizing single photon 
sources, which is of interest in present quantum informa- 
tion processing applications where the distinguishability 
and purity of single photons is of great importance. 
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FIG. 2: Simulated examples of the spectral density operator and corresponding measurement statistics, where Pa = Pa—Pb, for 
(top) a pure Gaussian distributed state, (center) a mixture of two Gaussian distributed states with different center frequencies, 
and, (bottom) a continuous mixture of Gaussian distributed states with different center frequencies (i.e. frequency jitter). 
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